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■ Abstract 

We consider complementary dynamical systems related to stationary 
Korteweg- de Vries hierarchy of equations. A general approach for find- 
ing elliptic solutions is given. The solutions are expressed in terms of 
Novikov polynomials in general quasi-periodic case. For periodic case 
these polynomials coincide with Hermite and Lame polynomials. As 
Q , byproduct we derive 2x2 matrix Lax representation for Rosochatius- 

' Wojciechowski, Rosochatius system, second flow of stationary nonlinear 

vector Schrodinger equations and complex Neumann systems. 

• y—t 

X 

H ■ 1 Introduction 

The problem of integrability and nonintegrability of Hamiltonian system with 
g degrees of freedom has been a subject of considerable interest for many years. 
Recently remarkable progress has been achieved in connection of study of sta- 
tionary soliton type equations. The method of restricted flows was introduced 
in [|l|, H 1^ as a non- linearization of the Korteweg-de Vries (KdV) spectral prob- 
lem and was generalized in 0| . The coupled Neumann system, the Neumann 
system, the Garnier type systems, the Rosochatius- Wojciechowski, Rosochatius 
and the Henon-Heiles type systems are examples of this type. In this relation 
important results about the algebro-geometrical interpretaion of these systems 
are obtained in B. General approach to completely integrable dynamical sys- 



tems of Neumann type is discussed in M, 10, h2|. Quasi-periodic solutions and 



spectral interpretation using both algebro-geometrical and spectral methods are 
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given in Quasi-periodic solutions (g — 2) and periodic solutions associated 
to Lame and Treibich-Verdier potentials are obtained in ^ . New method of 
constructing elliptic finite-gap solutions of the stationary KdV hierarchy, based 
on a theorem due to Picard, is proposed in |2^, |2^, |2^ . 

In this paper we are concerned with the following different approaches 

• integrable dynamical systems related to Hill's equation in the case of finite- 
gap potential |, i || | 0, 

• method of stationary flows and restricted flows 

• method of non-linearization of the KdV spectral problem ||^ 

• method of separation of variables |5^, ^ 

• algebro-geometrical construction [|[ |[ . 

We give unified construction based on algebro-geometrical approach. New 
solutions in terms of Novikov and Gelfand-Dickey (GD) polynomials are given 
in explicit form. 

The paper is organized as follows. In section 2 we construct the stationary 
flows associated to the KdV hierarchy strictly following |3^, ^ |4^. Some 
maps between completely integrable dynamical systems are presented. In sec- 
tion 3 we analyse the relation between restricted KdV flows and Garnier type 
systems following The map between stationary and restricted KdV flows 
are given in section 4 |42|. The list of known dynamical systems related to 
stationary KdV equations are presented in section 3. 

In section 4 we formulate Baker- Akhiezer function approach from algebro- 
geometrical and spectral point of view respectively. New solutions of integrable 
dynamical system associated to stationary KdV equations in terms of Novikov 
polynomials are presented . In particular these polynomials associated to Lame 
potentials coincide with Hermite polynomials and Lame polynomials. In section 
5 (2 X 2) Lax representations of Garnier, Rosochtius I and stationary second 



flow of vector nonlinear Schrodinger equation are given. Following pO] , 59 
r-matrix approach of these systems is discussed. 



2 KdV hierarchy and Gelfand-Dickey polynomi- 
als 

In this section we follow the geometrical construction of |3^, ^ ^ 
with only little changes of signs and different spectral parameter. Let M be a 
bi-Hamiltonian manifold: if the associated Poisson operator P'*' := Pi + AXPq 
admits as a Casimir a formal Laurent series h{X) 

h{X) := hjX-i (2.1) 

j>0 
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then ho is a Casimir of Pq and the coefficients hj {j > 0) are the Hamiltonian 
functions of a hierarchy of bi-Hamiltonian vector fields Xj : 

Xj = Pidhj = Podhj+i, (j > 0) (2.2) 

At any point u G M, the bi-Hamiltonian flows are given by ^ = Xj{u), tj being 
the evolution parameter of the jth flow. The vector field (2_^) are Hamilto nian 
also with respect to the Poisson operator P"^. In fact the recursion relation (2.2) 
can be written as 

Xj = P^d/i(j')(A), /i(^')(A) {\'h{\))+ (2.3) 

where the index + means the projection of a Laurent series onto the purely 
polynomial part. 

Let M be the algebra of polynomials in u, Ux, Uxx, ■ ■ ■ (u — u{x) is a C°° 
function of x and the subscript x means the derivative with respect to x), and 
le Pq and Pi be the two Poisson operators of the KdV hierarchy |45| : 

d d^ d 

Po--=-r^ Pi — - 4ii— - 2m,. (2.4) 

dx dx'^ dx 

The gradients of the Casimirs of the associated Poisson operator P'*' can be 
obtained searching the functions v{X) :— X]j>o '^i-^ 'which are solutions of the 
following equation: 

B\v{X),v{\))^aiX) (2.5) 

where a(A) = X]j>o '^J'^ ' ^-"^^ constant parameters and P'*' is the bilinear 
fimction 

B^{wi,W2) := WixxW2 + WiW2xx " WlxW2x " 4(m - X)wiW2. (2.6) 

In fact is related to the Poisson operator through the relation 

-^B^{W1,W2) =WiP^W2+W2P^Wi. (2.7) 



Equation (2.5) can be solved developing the left-hand side as a Laurent series 

P^«(A),t;(A)) = ^ PfcA-^ (2.8) 
fc>-i 

so that, for each a(A), it furnishes the cocfhcicnt of the solution v{X) (unique 
up to a sign). The solution corresponding to a(A) = —A is the so-called basis 
solution v{X); its first coefficients are 

vo = l, vi = l, v = 2{uxx + Su^), (2.9) 

V3 = 2{Uxxxx + 5ul + lOUxxU + lOu^), (2.10) 
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and so on, namely the gradients of the first KdV Hamiltonians. In what follows 
we shall consider also the function 7;(A) = c(A)w(A), which is a solution of ( ^.5[ ) 
for 

a(A) = -Ac2(A), c(A) = 1 + XI '^i-^^^ (^-H) 

where the coefficients Cj are free parameters. In this case the first 1-forms of 
the hierarchy are vq = l,vi ~ Vi + Ci,V2 = V2 + ciVi + C2, and so on. 

The coefficient Bj, can be expressed through the GD polynomials. For each 
Laurent series v (A) let us consider the functions B'-'^^X) B^{v{X),v'^''^Xj), 
where v'-^^X) :— (A'^w(A))^; these functions have the form 

fc-i 

5(fc)(^) = X'^+^I + J2 A'"^(Po, + vovj+i) + X X-'p,k. (2.12) 
i=i j>o 

It can be shown that 

B-i = -Vg, Bk = PQk - VQVk+i. (2-13) 



Furthermore, if v{X) is a solution of (2.5), the coefficients pjk in (2.12) are 
polynomials in u and its x-derivatives. They will be referred to as Gelfand- 
Dickey (GD) polynomials and the function B^ as their generating function. 

The fundamental property of the GD polynomials, stemming from (2.5), 
( p. 12 ) is the following relation with the gradients vj = dhj and the bi-Hamiltonian 
vector field Xk'. 

—pjk^VjXk. (2.14) 

We report some GD polynomials to be used in what follows {vq — 1): 

Poo = 4u - wi , 

Poi = 8uvi-vl-V2 + 2vixx, 

Pq2 = 4uUi + 8UV2 - 2viV2 -V3- vl^ + 2viVixx + 2v2xx, 

(2.15) 

p^2 = SuViV2 - V2 + "iuVs - V1V3 - V4 + 2vixV2x +2v2Vixx +2viV2xx + Vsxx, 

Pkk 2vkxxVk ~vl+ Auvl- 

The GD polynomials corresponding to the basis solution t'(A) are the polyno- 
mials defined in proposition 12.1.12. 

2.1 The method of stationary flows 

The method of stationary flows was developed in order to reduce the flows of 
the KdV hierarchy onto the set Alg of fixed points of the gth flow Xg of the 
hierarchy 



M, := 



{u\Xg{u,Ux,...,u^^a+i))^QY (2.16) 
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As Mg is odd-dimensional it can not be a symplectic manifold; nevertheless we 
will show that it is a bi-Hamiltonian manifold; it will be referred to as extended 



phase space. Moreover, Mg is naturally foliated, on account of (2.2) and (2.4), 



by a one-parameter family of 2g-dimensional submanifolds Sg given by 



Sg := |w|ug+i(u, Mx, ■ • • , M*^^-' = c| 



(2.17) 



(c being a constant parameter), which are invariant manifolds with respect to 
each vector field of the KdV hierarchy, due to the invariance of the functions v^. 
So Mg can be parametrized naturally by "yi, . . . , Vg+i and by their x-derivatives 
Vix, ■ ■ ■ , Vgx- We shall use these coordinates in what follows. 

From the computational point of view, one proceeds as follows. 



(i) Due to (2.3) and ( |2.5| ), the manifold Mg is defined by the solutions u of 
the equation 

B^{v{X),v^''>{X)) = X^aiX) (2.18) 



where v{X) = J2^j=o'^j'^ = S^l-i '^j'^ particular, if a(A) 

— Ac^(A), Mg is given by 



u\X„ 



7 . '^j-^n-j 
J = l 



(2.19) 



i.e. by the solutions of the Lax-Novikov equations. Taking into acount (2.12) 
and choosing a_i = —1, by equation (2.18) the coefficients of A^+^ we get 



Vq = 1; from now on we put vq — I- Moreover, equating the coefficients of the 
other powers of A we get the following system: 

PQk-Vk+i=ak, (/c = 0. . . ,n - 1) pj„ = a„+j, (j = 0, . . . ,7i). (2.20) 

(ii) In order to obtain the first Poisson tensor Pq: we eliminate u = vi/2 - 



ao/4 from ( ^.2C ) using the first equation (fc = 0) and we extract the system of 
n second-order ODEs in the Vj (j = 1, . 



,n 



Pok ~ Vk+i ^ ak, (fc = 1 . . . , g - 1) pog = ag 



(2.21) 



which will be referred as Po-system. The remaining equations (|^ will furnish 
a set of g independent integrals of motion. In order to obtain a second Poisson 
structure, we consider the following system: (Pi-system) 



Pofc - Wfc+l = flfc, (fc = 5-1) Pgg = a2g: 



(2.22) 



with u as above. 

(iii) The Pp-system (2.21) and the Pi-system (2.22) can be written as canon- 
ical Hamiltonian systems 



(0) 



rkx 



dsk 



Skx 



dH, 
dvk 



(0) 



(2.23) 



5 



Qkx = -TT—, Pkx = TT— (2-24) 

opk oqk 

have n integrals of motion given by 

K-j = -^Pjglr = flg+j, (j = 1, • • ■ , ff), (2.25) 

Hj = -^Pjglx = ag+j, (j 0, . . . ,g - 1), (2.26) 

Moreover, the map $ : Mg Mg in the extended phase space generates a 
second Poisson structure. 



3 Integrable dynamical systems related to hier- 
archy of stationary KdV equations 

In the recent years, remarkable progress has been achieved in the description 
of those quasi (periodic) potentials which belong to a given spectrum. Many 
integrable systems of differential equations are shown to be closely connected 
with Hill's equation in the case of a finite gap potential. The coupled Neumann 
system. The Neumann system, and the Rosochatius systems are examples of 
this type. 

In this paragraph we are concerned with the following completely integrable 
systems 

The Garnier system 

The g-dimensional anisotropic harmonic oscillator in radial quartic potential, is 
obtained when = rji i — 1, . . . , g 

^^xx = |^2^e|+a^j (3.2) 

Another interesting integrable system was proposed recently ji^ , we call it the 
Rosochatius- Wojciechowski system. In our context, this system is obtained by 
the Deift elimination procedure. Let 

S,, = tpi exp(6lj), = exp(-6'j), ^/f ^ = iplOi,^, (3.3) 



then equations (3.1) transform to the Rosochatius- Wojciechowski system 



i^^xx = I 2 ^ V'l I V'z + a^Vz - hl^l « = 1, • • • , 5 (3.4) 



6 



and the Hamiltonian is given by 
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^ n g 



^ = Ex?-(E^n -E«.^f-E4 (3.5) 

j=l \k=l I 3 = \ 3 = \ ^3 

where Xj = V'ja; are canonical momenta, 
(ii) The coupled Neumann system 



lixx + I 2 E ^i^ii^ii + iixf]rx ) = hli, (3.6) 
?7ix2: + I '^'^bj^jfij + £,jxV]x ) % = bifji. (3.7) 

with constraint X]f=o ^^'^^ = 1' where bo < bi < . . . < bg are fixed real numbers. 
The Neumann system is obtained when = rji 

Ixx + M E ^^4" + j = (3-8) 



This system describes the motion of uncoupled harmonic oscillators £^ixx = ^iCi 
constrained by the force X]?=o(^j^4^ +^i^x) move on the unit sphere X]i=o 



1. Let 

= ^jji exp(6ii), % = -(Ai exp(-^,,), \//- = i^'i^ix, (3.9) 



then by Deift procedure (|3.9|), equations (3.7) transform to Rosochatius system 

r~„. = - f E ^^^1 + ^1. - 4 - S + (3-10) 
where Ef=o ^? = 1- 



4 Baker- Akhiezer function 

We review in this section some basic facts about Baker- Akhieser function which 
will be used in the sequel. 

Let K be the hyperelliptic Riemann surface /i^ — n?f:o('^~'^«) ~ R{^)- The 
points of K are pairs P = (A, R) and \{P) is the value of the natural projection 
P — > A(P) of K to the complex projective line CP^. 

For given nonspecial divisor D, there is an unuque Baker- Akhiezer (B A) func- 
tion ^{t, A), such that 

(i) the divisor of the poles of ^' is Z?, 
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(ii) ^ is meromorphic on K\(x 

(iii) when P — > oo 



^'(a;,P)exp(-fcx) = 1 + ^ms(x)fc~^ (4.1) 



is holoniorphic and k — A(P) is a local parameter near P = oo. 
There is a unique function u{x) such that 

4-,, - = A(P)*, (4.2) 

where 5' is a BA function. Inserting expansion (0) into (U), we obtain 

- 2mi,(x)^ - A(P)* - exp(fca;)0(fc-i), (4.3) 



and due to the uniqueness of we prove (4.2), with u{x) = 2mix{x). 

By the Riemann-Roch theorem, there exists a unique differential and a 
nonspecial divisor D'^ of degree g such that the zeros of are D + D'^ and the 
expansion at P = oo, Vt{P) = (1 + 0{k~^ j)dk. 

For given nonspecial divisor D'^, there exists a unuque dual Baker- Akhiezer 
(BA) function such that 

(i) the divisor of the poles of ^' is D'^ , 

(ii) ^' is meromorphic on K\oo 

(iii) when P — > oo 

*^(a;, P) exp(-fca;) = 1 + XI W^"": (4-4) 

s=l 

Fix r to be the hypereUiptic involution P = (A, P) ^ (A, — P), then we have 
P)^ = tD, ^-^(2;, P) = ^'(x, rP). Let ^f^^ /i,(0) be the A-projection of D, and 
be the A-projection of the zero divisor of '^{x,P). The function 
\l/(t, P)^''^(t, P) is meromorphic on CP-*- and the foUowing identity takes place 

^{x,P)^^x,P)^^^ (4.5) 

where F{x, A) = ni=i(-^ ~ Introduce the Wronskian 

{*(a;,P),*^(a;,P)} = *^(a;, P)*^(a;, P) - *(a;, P)*;(x, P) = (4.6) 



2VP(A) 



nLi(A-M.(0)' 
and the differential il is given exphcitly by 

^{P) - ^ - M.(0))/v^rfA. (4.7) 



2 



8 



We assume that E{P) is a meromorphic function on K with g + I simple poles 
oo,pi, ■ ■ ■ ,Pg and at P — > oo, E{P) — fc + . . ., and E{P) is meromorphic function 
with 5 + 1 simple poles qo,qi, . . . ,qg and at P ^ cx), E{P) = k^^ + . . .. We also 
suppose that the divisors of poles of E{P) and E{P) are different from D, D'^ . 
Let 

^^fj^ = Resp^g^En, 6, = Afe), * = 0,...,5 (4.8) 
= Pesp=,.^^i, a, = A(k), z = (4.9) 



then 



\i=0 J i=0 

a 

(x) = 2^^iT]i + const. (4.10) 



Let us construct the meromorphic differential E'ii'i''^Cl. By direct computations 
we have 

a a 

Resp^g^ E^-^^n + Resp=^ E'i^^n = ^ l.fy, - 1 = 0, (4.11) 

i=0 1=0 

where = Resp=g. Eil. Differe ntiat ing Y^l=o ^i^i — 1 twice and using Eq. 
(^), we obtain first expression in ( 4.10| ). The eigenvalue equations 

^ + U{x)%, (4.12) 

flixx = (Hli) + u{x))f]i, (4-13) 



by replacing u{x) from (4.10) are the coupled Neumann system (3.7). By com- 
putations of the same kind, we have 

a ^1 
^ Resp=p. E'i'-ii^n + Resp=oc E'i'-^^n = ^ ^^r/i - u{x) + -const. = 0, 

1=1 j=i 

(4.14) 

where ^iij^ = Resp=p. EQ,. The corresponding eigenvalue equations are the 



Garnier system (3.1 



4.1 Spectral interpretation 

Let p be a positive real divisor of degree (/ on a real hyperelliptic curve 

^2 = i?(A) = [](A-A,), Ao <Ai <,...,< A2g. (4.15) 

i=0 
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The projection X{pi) lie in the closed lacunae [X2i-i, A2i]. The following consid- 
erations, due to Jacobi, can be used to construct such a divisor. 

Each divisor p determines and is determined by a system of polynomials 

i(A) = n(A-A(ft)), C'(A) = i(A)^ "{/^ 

B{X) = A^+i + . . . , (4.16) 

of degrees g, g — 1, g + 1, respectively, with R = — AB. The complementary 
divisor q is also determied by this construction. This is the content of the 
following step. 

For a given spectral data 

Ao = 0<Ai <,...,< Aag, A(p^) e [A2i-i,A2i], i = 1,...,5 (4.17) 
thete exists 

A(gj), z = 0, ...,.g A((7o) e (-oo,Ao], \{qi) e [\2z-i,\2i\, (4.18) 

such that R = — AB, and the projections A(ft) are the roots of B. 

Note that the functions E{P), E{P) are meromorphic on K and the following 
formulas are immediate 

E{P) = (/R(A) + (7(A))/i(A), C{pi) = ^/m^ (4.19) 



E{P) = (yi?(A) + C'(A))/S(A), C{q,) = -/rR (4.20) 

Now we recall some facts from the periodic theory of Hill's equation. We 
suppose that u{x) is a real finite-gap potential, i.e. the operator L has only 
2g + I simple eigenvalues Ao < Ai <,...,< A2g and the rest of the spectrum 
consists of double eigenvalues. The periodic spectra of L is determined by the 
combined eigenvalues of the periodic 

L^2^= \2^h^. f {x + 1) = f {x) , i = 0, . . . , g (4.21) 

and the antiperiodic 

L f2^-l = X2^-l .f2^^1, f {x + I) = - f (x) , I = 1, . . . , g (4.22) 

eigenvalue equations. The intervals (— oo, Aq], [A2i-i, A2i] are termed lacunae. 
The Floquet solutions(periodic BA function) and the corresponding Floquet 
multipliers, are given by 

^(x, A) = [F{x, A)/F(0, A)]^/^ exp (^J^ .JR{X)/F{x', X)dx'^ , (4.23) 
*(x + l,A)=p+(A)vE'(x,A), (4.24) 
*^(x,A) = [i^(x,A)/F(0,A)]^/^ exp(^- ^J R{X)/ F{x' ,X)dx'^{^.2b) 
^^{x + l,X)=p-{XW{x,X), p± = exp(±p(A)), (4.26) 
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where 



p{X) = / y'R{X)/F{x, X)d2 



(4.27) 



Note that if A is in the periodic spectrum, ^'(x, A2i) — f2i, i — 0, . . . , g is a 
periodic eigenfunction, and ^'(a;,A2i-i) = /2j-i,* — ^,---,g is an antiperiodic 
eigenfunction. It is well known that the projections of the zeros of the Floquet 
solution define the auxiliary spectrum of L. 
The following expressions hold 



(4.28) 



a 



n(A(K)-A^,(0))/i'(p.), 



Un^ = \{{Kq^)-Ni^))|B\qi), 



(4.29) 



a 



J{{\{q,)-^iM)|B'{q^)■ 



Using ([f.iq),([4.20D we obtain 

Resp=p^E{P)'^'^^VL 



where ^°77° is given by (4.28). 
Let 



-2i 



nf^,=o (A24-M,(o)) 



1 /2i — 



f{{X{p,)-^,,{x))/A'{p,), 



IIU (^2.-M,(t)) 



m^,=o (A2,-M,(o)) 



= *"(A20, (4.30) 



i — . . . , g, then the expressions ( 4.101 ), (4.29) are the famous McKean-Moerbeke 
expansion of the potential u{x) in terms of squares of the eigenfunctions 

(g a g \ 

2^1 A2^/2Ve24 + ^24-1 - ^2.: + Ao , (4.31) 
1=0 1=1 i=l / 

where the following identity among the squares of eigenfunctions hold on 

g 



(4.32) 



4=0 



The results of this section, may be summarized by following: 
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Let u{x) be a real nonsingular finite-gap potential. There exists g eigcn- 
functions \E'(y»i), . . . , "^{Pg) and g + I eigenfunctions "^{qo), • • • , of Hill's 

equation, corresponding to the eigenvalues A(pi), . . . , X{pg) and A(go)j ■ ■ ■ , A(gg), 
respectively, such that 

(i) 

g g "^g 

u{x)=2Y, M'(pO*"(Pi)e-' + 2^A(pO -E^^' (4.33) 

i=l i=l i=0 

g 

(ii) 

9 9 "^9 

u{x) = 2 £ A(ft)*fe)*"(9i)er' - 2 £ A(<j,) + £ A^, (4.34) 

i=0 i=0 i=0 

g 

= ll(^{Qi)-Nm/B', ^iqi) = ^{t,X)\x=,„ i = 0,...,g, 

g 

^e-2*(gi)*-(ft) = l. (4.35) 

i=0 

The corresponding eigenvalue equations are the Garnier and coupled Neumann 
system. 
Let 

g g 

4-1 = l[i^2i-i - X2j-i)/ n(^2i-i - nm, 

g g 
fl-i = ]Ji>^2i-i - Hj{x))/Y[{X2i-i - fij{0)), i = l,...,g 

then we have the following expansion of the potential ^(a;) in terms of squares 
of antiperiodic eigenfunctions 

g g 2g 

u{x) = 2^ fl-ie^l, + 2^ A2.-1 - J2 (4-36) 

i=l i=l i=0 

We call the dynamical systems such in (i) , (ii), complementary dynamical 
systems. 
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4.2 Solutions in terms of auxiliary spectrum of Hill's equa- 
tion 

The solutions of the Gamier system in terms of auxihary spectrum fJ.j{x), j = 
1, ... ,5 are 

= ^°[F{x, a^)/F{0, a,)]i/2 exp (^J^ ^R{a^)/F{x', a,)dx'^ , (4.37) 
7], = Tj^[F{x,a^)/F{0,a,]^/^ exp ^R{a,)/F{x' ,a,)dx'^ ,{4.38) 



where fj,j (x) satisfies the following system of differential equations 

A ' ' 

dx 



— /ij(a;) = 2^R{iij)/Y[i^j.jix) - fik{x)), 



with initial conditions 

g 

M,(0) G [A2,-i,A2.], e°r?° = F(0,a,)/nK 
Differentiating expressions 



and 

we have 



T(x,P)U=Afe) = ^loge»(i) 



,1 d 



(4.39) 
(4.40) 

(4.41) 
(4.42) 



(4.43) 



,1 



i=i 



(4.44) 



direct integration of ( [4.43| ) , ( |4.44| ) gives the solution s ([4.37] ) , ( [4.38[ ) . The function 
T(x, P) has (7 poles at /Xj(x), then the numerator of (4.43) is zero when A = l^j{x) 
and the following system takes place 



(4.45) 



13 



This is another form of the system (4.39). In the same way, we can obtain the 
solutions of the coupled Neumann system by replacing X{pi), i — 1, . . . , g with 
A(g,), * = 0,...,5in (Q, ( p^ ). 

Let X{pi) be the antiperiodic eigenvalues A2i-i, i — l,...,g. Then the 
exponential function in (4.37), (4.3S) cancel, = rfl and the solutions of the 
(^-dimensional oscillator are 



- n (A2.-1 -M,(x))/n(A2.-i - \2,-i) 



(4.46) 



Let X{pi) — Ui be in a general position, i.e. X{pi) G [A2i-i, A2i] and by Deift 
elimination procedure we may identify with $^f[F{x, ai)/F(0, Oi)]^/^ and 

e^ = r /RK)/n(a. - ^i,{x'))dx' , = (4.47) 

and, hence, the solutions of the Rosochatius-Wojciechowski system are 



(4.48) 



Inserting explicit expression of BA- function given by ( I.23| ) in Hill's equation 
we have 

^F,Ax,X)F{x,X) - \f^{x,X) ~ {u{x)+X)F\x,X) = -R{X). (4.49) 

The polynomial solution of ( [1.49 ) below we will call Novikov polynomial |p5| . 
Assuming that Novikov polynomial dependts on time t, the zero curvature rep- 
resentation for KdV hierarchy of equations have the following form 



Mt(A')-L.(A') + [Af(A'),L(A')], 
where matrices L and M are given by 

L(A') = 



(4.50) 



M'(A' 



-F,(x,t,A')/2 F{x,t,X') 

-FrAx, t, A')/2 + Qix, t, X')Fix, t, A') F, [x, t, X')/2 

1 

Q{x,t,X') 



The equation ( 4.50 ) is equivalent to 



dQ 
dt 



= -2 



-^dl~Q[x,t,X')d^~]^QAx,t,X') 



Fix,X'). 



(4.51) 



where Q{x,t, A') = u{x) + A' in the case of KdV hierarchy. Equation (4.51) is 
called the generating equation. For a different choices of the form of F{x,t, A') 
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and Q{x,t,X'), this procedure leads to different hierarctiies of integrable equa- 
tions, as an example to the KdV, nonlinear Shrodinger and sine-Gordon hierar- 
chies or to the Dym hierarchy |Q . The Lax representation = [M, L] yields 
the hyperelliptic curve K = (//', A') 

Det(L(A') - Ai'I) = 0, (4.52) 

= -Iff,, + + (A' + u)F = R{X'), 

generating the integrals of motion for stationary KdV hierarchy. 



Using the equation (4.49) and the following expansion of potential u{x) in 
terms of squares of eigenvalue functions £,'^{x) = l3k{x) 

9 9 '^9 

u{x)^2Y,ekix) + 2Y,'^^-Y.^^^ (4.53) 

i=l i=l k=0 

have the form 

lPk.A - ]PI - {u{x) + ak)l3l = = -h (4.54) 

where we use the solutions F{x,ak)/akj of Rosochatius-Wojciechowski system 
and flfcj = Yll-ij{ak - aj). Denoting fk = R{ak)/akj and d = X^Li ~ 
i X]fe=o -^fc ^'^^ original variable /3j. — we have the following equation 




^k..^2{yjl+d]^k-^. (4.55) 



To understand the role of GD polynomials and of their generating function 
in the construction of a map between stationary and restricted flows of KdV 
equation and exact solution of completely integrable systems related to Hill's 
equation let us consider the following system: 

Poo -vi= ao, Po ^1 - ^ , (4.56) 

P^''f3k = 0. (fc = l,...,n) (4.57) 

whereAi, . . . , A„ are fixed parameters, P^'' := Pi + 4Aa;Po (-Po and Pi being 
the two KdV Poisson operators). This is a system of {g + 2) equations in 
u,vi,(3i, . . . ,f3g. The second equation will be referred to as the Po-restriction 
of the first KdV fiow Xq — PqVi = vi,, and the last n equations define the 
kernel of g Poisson operators extracted from the Poisson operator. On account 



of (2.15), (2.4) and (2.6) this system is equivalent to the following one: 



y + T' «i=E/5j+c> S'H/3fe,/3fc) = /fe, (4.58) 
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where c and fk are free parameters and is just the generating function of 
the GD polynomials. 

Using the first two equations to eliminate u and vi from the last g equations, 
one gets a system of n second-order ODEs for /3i , . . . , /3g : 

2f3k.A - PL + 2Pl + - ^^^k = fk, k = l,...,g (4.59) 

where d :— c + ao/2. Introducing the so-called eigenfunction variables ip'j = 
Pj and the momenta Xj — ^jxi equations (4.59) can be written in canonical 
Hamiltonian form 

V3^^^ — ' Xj^ = -^— , J = l,...,n (4.60) 

with Hamiltonian 

n / n \ n n f. 

^« = Ex|- E^n -E«.^I-E4 (4-61) 

3 = 1 \k=l I j=l j=l ^3 



A set of integrals of motion is 



E- 



k=l / ^3 
1^1 fk^^ 



k^3 



4 + 

+ {i^jXk - i'kX3?\ ■ (4-62) 



where we denote ajk = aj — Ofc. 

Now we shall construct a map between the g-th stationary flow and the previ- 
ous restricted flow of the KdV hierarchy. To this end we extend the correspond- 
ing phase spaces, regarding some free parameters in the Hamiltonian functions as 
addiitional dynamical variables. As for the Pi -formulation of the stationary flow 
we extend its phase space to a {<ig+ l)-dimensional space, Af„, with coordinates 
(ftjPfe; flO: ■ ■ • 1 Og-i, o,2g)', aualogously we consider the Po-formulation of the first 
restricted flow in the extended space Mg with coordinates {ipk ,Xk', fi, ■ ■ ■ , fk,d). 

Let us consider the solutions qk of the dynamical equations (2.24); then 
ti(£')(A) given by 

i;(^)(A) = A(g2(A))(s-i) (4.63) 



^9 



with q{X) = 1-1- X]f=i Ij-^ satisfies (2.17), and consequently satisfies the 
following equation: 

B^(t;(f), «(")) = A29d(A), (4.64) 

where, as above, we put u — i;i/2+ao/4. So, for each g-tnp\e of distinct complex 
parameters aj, any solution ^^^•'(A) fulfills the system 

B''^{v^3\ak),v^3\ak)=al'd{ak), {k = l,...,g) (4.65) 
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where w^^-'(afc) :— v^^^x^a^- In order to have a solution also satisfying the 
second equation vi — X)f=i Pj + c, the Lagrange interpolation formula can be 
used. It allows us to represent the polynomial u"(A) by 

^(")(A)=a(A) (l + ^^^ j , (4.66) 



where a(A) = nj=i(-^ ^ "j)' ^^"^ 

f3j-^^^^^,ik^l,---,9)- (4.67) 

(a'(A) means the derivative of a(A) with respect to A. 

Obviously the g functions Pk are solutions of the following system 

2pk..Pk - PL + 2Pl {j2Pj+d^- XkPl = fc = 1, . . . , .9 (4.68) 

Furthermore, Pk satisfy the so-called Bargmann constraint 

g 

- aj) = vi, (4.69) 



as one can verify by means of (4.66). 

The function S'^ is also a generating function of integrals of motio n for 
Garnier system. Indeed evaluating the function by means of ( 4.66) an d 
eliminating the first x-derivatives of Xk by means of Hamilton equations ( 2.24 ), 
one gets 

where Ij are the functions. Taking in this equation the residues a t A = Uj it 
follows that the functions Ij are integrals of motion along the flow ( 2.24 ). 
Let X{qi) — hi be in a general position, i.e. 

\{qi) = bi G (— oo, Aq], [A2i-i, A2i] and by Deift elimination procedure we 
may identify with |^[F(t, 6,)/i^(0, 6,)]^^^ and 

l-x 9 

~e.= ^/Rlbi)l\{{h- ^i,{x'))dx', (4.71) 
and, hence, the solutions of the Rosochatius system are 

g g 

= n (&- ^^M)/tl^b, - bj). (4.72) 

j=0 i^j 
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where i,j — 0, . . .g. Now wc illustrate the general aproach with some simple 
examples 

Example 1 Let u{x) = Qp{x + lo') be the two-gap Lame potential with 
simple periodic spectrum (see for example [p^) 

= -^3^2, Ai = -3eo , A2 = -3ei , A3 = -3e2 , A4 = \/ig2- (4.73) 
and the corresponding Hermite polynomial have the form 

F{p{x + Lo'),X) = ~ 3p(a; + uj')\ + 9p^{x + lo') - -02. (4.74) 
Consider the following genus 2 nonlinear anisotropic oscillator with Hamiltonian 

H - \{p\+pI) + \{ql + qlf - \{aiql + a^ql), (4.75) 

where {qi,pi), i — 1,2 are canonical variables with pi ~ qi^ and 01,02 are 
arbitrary constants. The simple solutions of these system are given in terms of 
Hermite polynomial 

.? = 2^, .2^ = 2^^, (4.76) 
A2 — Ai Ai — A2 

Let us list the corresponding solutions 

• Periodic solutions expressed in terms of single Jacobian elliptic functions 

The nonlinear anisotropic oscillator admit the following solutions: 

qi ~ Cisn{ax,k), (4.77) 
q2 C2cn(aa;,A:), (4.78) 

where amplitudes Ci ,6*2 and temporal pulsewidth 1/a of are defined by pa- 
rameters fli and 02 as 

cP'k^ = 02 — 01 

2 _ „ , „,2 „,2;2 



C( = 02+a -a fc^, 

r2 

2 



Cn = oi + a^, 



where < fc < 1. 

Following our spectral method it is clear, that the solutions 4.78| are associ- 
ated with eigenvalues A2 = —62 and A3 = —63 of one - gap Lame potential. 

• Periodic solutions expressed in terms of products of Jacobian elliptic func- 
tions 
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qi = Cdn{ax,k)sn{ax,k), (4.79) 
92 — Cd'a{ax,k)cii{ax,k), (4.80) 

where sn,cn, dn are the standard Jacobian elhptic functions | |5^ , k is the mod- 
ulus of the elhptic functions < /c < 1, an the wave characteristic parameters: 
amplitude C, temporal pulsewidth 1/a and k are related to the physical param- 
eters and, k through the following dispersion relations 



a" 



2(4q2 - ai) 
5 

(4a2 - ai) 
15 ' 
5(a2 - ai) 



4a2 - ai 

We have found the following solutions of the nonlinear oscillator 

qi — C k^ cn{ax T k)sn{ax , k) (4-81) 

92 = Ca^dn{ax,k) + Ci (4.82) 

where C, Ci, a and k are expressed through parameters ai and 02 by the 
following relations 
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a2 — fli 
C(4ai - 02) 
5 



2 



\Jj{a\- a\) + 002 — 3ai 



= ^(2a2-3ai-HY/^(a^-a2)). (4.83) 
Periodic solutions associated with the two-gap Treibich-Verdier potentials 



Below we construct the two periodic solutions associated with the Treibich- 
Verdier potential. Let us consider the potential 

u{x) = &p{x + uj') + 2 (£i_M(£L^ (4.84) 

p[X + LU'j - d 

and construct the solution in terms of Lame polynomials associated with the 
eigenvalues Ai, A2, Ai > A2 



Ai = 62 + 2ei + 2 V(ei - e2)(7ei + 2e2), 

A2 - 63 + 2ei + 2^(61 - e3)(7ei + 263). (4.85) 
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The finite and real solutions qi , q2 have the form 

qi — Cisn(z, A:)dn(z, fc) + C2sd(z, A:), 

92 = C3Cn(z,A:)dn(z,fc) + C4cd(z,A:), (4.86) 

where Ci, i = 1, . . .4 are constants and have important geometrical interpreta- 
tion The concrete expressions in terms of fc, Ai, A2 are given in [0 

Analogously we can find the elliptic solution associated with the eigenvalues 

Ai = 62 + 2ei + 2V(ei - e2)(7ei + 2e2), A2 = -6ei, (4.87) 

We have 

qi = CiAn{z,k) (4.88) 

52 = Cisn(z,A:)dn(z,fc) + C2sd(z,A:), (4.89) 

where C is given in . 

The general formula for elliptic solutions of genus 2 nonlinear anisotropic 
oscillator is given in 



N 



- ^ - ( 2A? + 2Ai p{x - 
A2 — Ai \ 



Ng2 



6 p{x-x,)p{x-Xj)-^+ J2 A,A, |,(4.90) 

i<i<j<N l<i<i<5 
1 / 

^ 2A^ + 2A2 Y - 



Ai — A2 



i<i<j<f^ l<*<i<5 / 

where Xi are solutions of equations X^i^^j p'(2;i — a;^) = 0, j = 1, . . . , and N is 
positive integer. 

Example 2 Garnier type system. Consider the following genus 2 Garnier 
type system with Hamiltonian 



1/ 2 I 2\ I -'^i' 2 I 2\2 1/ 2 I 2\ I /l I /z 

oil'i +P2) + 7(91 +92) - o(7i9i +7292) + -2 + 72' 
z 4 z yi 92 



where {qi,pi), i = 1,2 are canonical variables with = g^j^ and 01,02 are 
arbitrary constants. The simple solutions of these system are given in terms of 
Hermite polynomial in the following form (4.48) 



2_ i^(x,ai) ^_^F{x,a2) 



q^^^2^-^^^, (4.93) 
ai — 02 02 — oi 
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the same settings of periodic spectra, Lame potential and Hermite polynomi- 
als ( }4.74 ) as in example 2. The main results from the general theory are the 
following: 

i) the parameters /i and /2 are expressed in terms of algebraic curve by 
( p4| ) 

. _ R{ai) , R{a2) 

Jl — , J2 



v2 



R{X) = (A" - 352)(A + 3ei)(A + 3e2)(A + 863). 
ii) the parameters oi and 02 must lie in one or other of the intervals 

[-V33^,3ei], [362,363], [^35^.^) (4-94) 



and 7i = 3ai + 2a2, 72 = 2ai + 802. 

These results are in complete agreement with solutions obtained by different 
method in recent paper Q. 

Example 3 Simple solutions of the Henon-Heiles type system. 

We consider a generalized Henon-Heiles type system with two-degrees of 
freedom. 

Its Hamiltonian is 

Ho = \{pl +pI) + cil + Iq.ql + ^^+^(^ql + Iq^^ - ^q„ (4.95) 

where qi, q2,Pi,P2 are the canonical coordinates and momenta and 09,01,04 
are free constant parameters. This Hamiltonian encompasses the two cases 
Oq = 04 = and oo = oi = introduced in [ pO| . Moreover Hq is related with 
the Hamiltonian 

Hh = \{pI+ pI) + 9? + ^9192 + ^ + ^ (^9? + Bql) , (4.96) 
through the map 

q^=q^ + — -2B, 92 = ^2, oq = -2A + 125, ai = -A^ + 16 AB - 48B^ . 

(4.97) 

The function Hh is the Hamiltonian of a classical integrable Henon-Heiles sys- 
tem with the additional term 04/8^!. 

The function (4.95) is the Hamiltonian of the vector field obtained reducing 
X(){u) = Ux to the stationary manifold M2 given by the fixed points of the flow 
X2 + ciXi + C2X0 

M2 = -1- lOUxxxU + '^OUxxUx + iOUxU^ + Ci{Uxxx + 6UxU) + C2Ux = 

(4.98) 

where ci ~ — oo/2, C2 — — ai/2 -I- Oq/4. 

It can be obtained specializing to the case g = 2 the Hamiltonian of the 
Pi-system. In this case i?2^'' = Hq and the canonical coordinates and momenta 
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are, respectively, qi = ui/2, g| = — U2, Pi = lix, P2 = <l2x- The integrals of 
motion obtained by the reduction of the GD polynomials are 

-ffo = -\p02\x (4.99) 

H2^-\p22\x = ~'^ (4.100) 

Hi = -^Pi2U - (4.101) 

where -lpi2\x is given by 



P2qi -piP2q2 - 77'7i92 - + -7Z2 r9i92 + —92- (4.102) 



1 2 2 1 4 , 04'3'1 ao 2 , 2 

Next we will derive (2 x 2) matrix Lax representation for generalized Henon- 



Heiles system (4.95). Using Lax representation = [M, L], particular case of 
eq. ( 4.50| ) i.e. when there is no time dependence, we have 



F{x, A) = + IqiA - ^ql V = -F,/2 = -ipiA + ^q2P2, 
2 lb 4 lb 

W = -F,J2 + OF = A^ - i^qi + iao)A2 + 



^i'^i + 16'^^ ^ 16"^ + 8"°'^^)^ + 16^2 + 64^' ^) = ^ - «i - 



The corresponding algebraic curve have the form 

= A^ - ifloA* - — aiA^ + 8Ho\^ + 32 i?iA 5— 04. (4.103) 

^ 4 16 1024 ^ ' 



Using explicit expression for Hermite polinomial (4.74) we obtain the following 



simple solutions for the system (4.95) 



qi^-6p{x + Lo'), 92 = -16(9p(x + u;')'-^32). (4.104) 
where ao = 0, ai = 3.4.7(72, ^4 = — 3^.4'*g253. 

5 2x2 Lax representation and r-matrix approach 

The Lax equation for completely integrable systems discussed in the previous 
section 

L,(A') = [M(A'),L(A')] , (5.1) 
with matrices L and M given by 

r(y^ _ f V{x,X') U{x,\') \ 

- [ wix,y) -v{x,x') ) ^^■^> 



22 



is equivalent to the Garnier system, where U{x, A'), V{x, A'), W{x, A'), Q{x, A') 
have the form 
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Uix, A') = a(A') ( 1 - ^ ) , Vix, A') = -^U^ix, A')(5.4) 



A' - a,- / ' ' ' ' 2 



PF(x,A') = a(A') (A' + ^6r;. + ^|^] , Q(a;, A') = A' + 2 ^ e.'7..(5.5) 

Finally we point out one usefuU expression, which is easy to derive from Lax 
representation ( ^.iD 

W{x, A') = Uix, X')Q{x, A') - ^U{x, A')... (5.6) 

The Lax representation yields the hyperelliptic curve K = (/i'. A') 

Det(L(A') - //I) = 0, (5.7) 
generating the integrals of motion H, F^^\ i — 1, . . . , g. We have 

V^{x, A') + U{x, X')W{x, A'), (5.8) 



From (5.8) and explicit expressions of U{x, A'), V{x, A'), W{x, A') we obtain 
where 



_l_ {£,i^kx - S,kiix){riirikx - VkVix) 

Hi = ^ix^x - aiS,ir]i - ^^rji (^Yl^^kVkj , Ji = S.ix'm - ^iVix, (5.10) 

and X]f=i Hi is the Hamiltonian for Garnier system. Simple reduction rji = 
gives us the second flow of stationary vector nonlinear Schrodinger equation, 
where by * we denote complex conjugation. The complementary to the last 
system is complex Neumann system (see for example |]6^ ) 

i^xx + 2 6, ||,f + j 6: - (5.11) 

with J2i=o I'^jP — 1- Using the Deift elimination procedure we obtain new 2x2 
Lax pair for Rosochatius-Wojciechowski system. Bellow we list only the final 
results for Lax pair elements of considered in this paper dynamical systems: 
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• Rosochatius-Wojciechowski system 



U{x, A') = a(A') ' ^'^ = 



i=l 



A') = a(A') (a' + ^ V? + E ^^(V'ix - 4) ) , (5.12) 

\ i=l i=l ^ Ti J 

0(x,A')=A' + 2^^2_ 



i=l 

• Rosochatius system 



C/(x,A') = 6(A') ' ^(^'^') = 

A') = b{X') E ^^(C - ^) j , (5.13) 

9 

Q(x,A') =A' + 2EV'', 



i=l 

where 6(A') = n?=o(A' - ^i)- 
• second stationary flow of vector NLSE 



U{x, A') = a(A') - E j , y{x, A') = -^f/x(a;, A') 

W{x, A') = a(A') [a' + + , (5.14) 

\ i=l i=l '/ 

9 

2 



Q(x,A')=A' + 2El€i| 



• complex Neumann system 



u{x,\') = h{x') (e aTI^) ' ^^^'^'^ = -\u.{x,\') 

W{x,\')=a{\'){l-j^^\ (5.15) 



9 



Q(x,A') =A' + 2Ek.| 



i=l 
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Finally we want to point out that li for Rosochatius-Wojciechowski sys- 
tem coincide with expression given in (4.62). Another Lax equation have the 
following form 

L,(A') = [Af(A'),L(A')] , 
where matrices L and M are given by 

-F,{x,\')l2 F{x,X') 



(5.16) 



^^^'^ " * -F,,(a;,A')/2 F,(a;,A')/2 

V{x,X') U{x,\') 
W'{x,X') -F(a;,A') 



M'(A') - ( J ) • (5.17) 

where we made the following identification U{x,X') = F{x,X'). The Poisson 
bracket relations for the matrix L(X') |5^, ^ are closed into the following r- 
matrix algebra 

{Li(A'), W)} = MX', Ai'), ii(A')] - hi (A', m'), i2(/i') ] • (5.18) 
Here the standard notations are introduced: 

Li(A') -i(A')®/, L^ifi') = I ^ Lifi') , 

(5.19) 

ri2(A,M) = T^ r2i(A,M) =nri2(M,A)n, (5.20) 
A — /i 

and n is the permutation operator of auxiliary spaces. 

The Poisson bracket relations for the Lax matrix L{X') are preassigned by 
the initial symplectic structure. It is necessary to calculate only two brackets 

{Fix,X'),Fix,^^')}^0, (5.21) 

and 

{V{x,X'),F{x,fi')} = 



-F(x,A')F(x,m')E7T7T 



(A' -/ij(a;))(/i' - ^J.,{x)) 
F{x,X')F{x,fi') ^ / gjj gjj 



X' - fi' j-^yx'-nj{x) fx' - fi'^ix) J 

^ j[F{x,fi')-F{x,X')] , (5.22) 



A'-^i' 
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where we used a standard decomposition of rational function 

n 



and the following definitions 



g.7 



Res|,,^^^(,)F(x,A')-^(A) 



R-eS|;^=^^.(^)F ^(a;,A) 



(5.23) 



P]{x) = V{x, A)|A=^^.(,) = ^jR{fij{x)). (5.24) 

Another Poisson brackets may be directly derived from these brackets and 
by definition of the entries of the Lax matrix L{X) via derivative of the single 
fimction F{x, A') 



{V{x,X'),V{x,fi')} 
{W{x,X'),F{x,^l')} 







dx 



{V{x,X'),F{x,fi')} 



X' -fi' 



-[V{x,X')-Vix,^^')] 



(5.25) 



{W{x,X'),F{x,^i')} = 



{W{x,X'),W{x,fi')} 



-^^{n^,A'),F(x,M')} 
1 



A' 
1^ 

^2d^ 



-[W{x,X')-W{x,fi')] 



{V{x,X'),F{x,^l')} = 0. 



Applying the following transformation directly to the Lax representation 
L(A') we obtain a family of the new Lax pairs 



L'(A') = L(A') - a- ■ [</>(x, X')F{x, XT'] ^ , 
Af'(A') = M - ■ [(t>{x, X')F{x, X'y^^ - 





1 



(5.26) 



1 

Q{x,X') 



Here <j){x,X') is a function on spectral parameter and [zj^v means restriction 
of z onto the adjj-invariant Poisson subspace of the initial r-bracket. For the 
rational r-matrix we can use the linear combinations of the following Taylor 
projections 



[z]n -- 

or the Laurent projections. 



+ 00 



N 



(5.27) 



N 



k=0 
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New Lax matrix L'{\') js^, ^ obeys the linear r-bracket, where constant 
rij-matrices substituted by r^^ -matrices depending on dynamical variables. 

fi2(A',M') = ri2 



(A' - fi') 



(5.28) 



6 Conclusions 



In this paper we have given new exact solutions of the physically significant 
completely integrable dynamical systems. These solutions can be interpreted as 
eigenfunctions of suitable differential operator. New example of complementary 
dynamical systems (stationary second flow of the vector nonlinear Schrodinger 
equation and complex Neumann system) is presented. 
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